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Abstract 

A  turbine  stator  profile  that  is  analytically  predicted 
to  exhibit  good  aerodynamic  performance  over  a  wide  range  of 
transonic  operating  conditions  has  been  designed  for  the  Air 
Force  Aero  Propulsion  Laboratory’s  Heat  Transfer  Facility. 
Analytical  Turbine  Design  System  computer  programs  developed 
by  General  Electric  Company  were  used  to  generate  blade  pro¬ 
files/coordinates  and  perform  2-Dimensional  cascade  through- 
flow  calculations  (streamline  curvature  technique) .  Two 
aerodynamic  performance  standards  were  applied  to  pressure 
distribution  predictions  of  the  blade  contours  produced  from 
a  parametric  variation  study;  1)  equal  flow  expansion  over 
both  the  suction^ and  pressure  surfaces,  and  2)  minimum  suction 
surface /diffusion.  The  influence  of  the  following  geometric 
parameters  upon  the  design  of  a  transonic  turbine  stator  were 
investigated;  axial  location  of  maximum  section  thickness, 
stagger  angle,  leading  edge  bluntness  factor,  and  trailing 
edge  taper  factor.  Results  of  this  study  show  that  a  good 
transonic  turbine  stator  design  must  have  a  convergent- 
divergent  flow  passage  and  a  flat  suction  surface  aft  of  the 
throat  location. 


THE  DESIGN  OF  A  TRANSONIC 


GAS  TURBINE  STATOR 


I.  Introduction 


Background 

Modern  high  pressure  gas  turbines  operate  in  temperature 
environments  in  excess  of  1055°K  (1900°R) .  Due  to  the  thermal 
limitations  of  existing  materials,  the  requirement  for  blade 
cooling  has  become  an  essential  feature  of  gas  turbine  design. 

The  mission  objective  of  the  Air  Force  Aero  Propulsion 
Laboratory's  Heat  Transfer  Facility  (Figure  1)  is  to  test  and 
evaluate  high  temperature  turbine  blade  cooling  schemes  for 
advanced  turbopropulsion  applications.  The  current  test 
section,  a  rectilinear  turbine  stator  cascade  (Figure  2) , 
was  designed  during  the  mid-1960's  using  a  subsonic  turbine 
blade  design  philosophy.  Subsonic  blade  shapes  are  charac¬ 
terized  by  high  aspect  (span/chord)  ratios,  and  relatively 
thick  profile  shapes. 

Present  day  gas  turbines  are  designed  to  function  at 
pressure  ratios  which  cause  the  working  gas  medium  to  ap¬ 
proach  sonic  flow  conditions.  These  designs  are  motivated 
by  the  high  overall  thermodynamic  efficiency  possible  with 
a  transonic  pressure  ratio.  Small  cross-sectional  areas  and 
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Figure  1  Heat  Transfer  Test  Facility  Schematic 


Figure  2  Schematic  of  Test  Section  Cascade 
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thin  trailing  edges  are  generally  required  to  achieve  accep¬ 
table  aerodynamic  performance  in  a  transonic  flowfield. 

These  requirements  must  be  balanced  with  the  need  to  permit 
interior  cooling  systems  while  allowing  the  blade  to  retain 
mechanical  strength  at  high  temperatures. 

Because  the  heat  transfer  loads  and  thermal  stresses 
associated  with  a  transonic  profile  are  substantially  more 
severe  than  those  associated  with  a  subsonic  blade  shape, 
relevant  test  conditions  can  no  longer  be  established  in 
the  AFAPL  Heat  Transfer  Facility.  The  facility's  mission 
objective  is  seriously  compromised  until  a  new  transonic 
turbine  stator  contour  that  simulates  modern  gas  turbine 
environments  is  designed. 

Problem 

The  basic  intent  of  this  study  is  to  design  a  modern 
turbine  blade  shape  that  exhibits  a  high  level  of  aero¬ 
dynamic  performance  characteristics  in  a  transonic  flow 
regime.  Flow  through  a  transonic  cascade  passage  is  subsonic 
at  the  cascade  inlet,  choked  at  the  flow  passage  throat 
(minimum  area  location) ,  and  slightly  supersonic  (M  =  1.17) 
at  the  exit  plane.  Two  tasks  associated  with  the  accomplish¬ 
ment  of  this  design  objective  are:  1)  the  establishment  of 
aerodynamic  performance  criteria  based  upon  inviscid,  shock- 
free  flowfield  calculations,  and  2)  the  development  of  an 
analytical  design  optimization  procedure. 
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The  major  analyatic  design  tool  that  will  be  used  to 
generate  and  evaluate  various  balde  contours  is  a  set  of 
computer  programs  entitled,  "Turbine  Design  System”  (TDS) . 
Because  the  codes  were  recently  developed  and  have  not  been 
sufficiently  verified  for  accuracy  and  consistency,  a 
secondary  problem  of  verifying  the  TDS  programs  arises.  The 
computer  codes  were  generated  by  General  Electric  Company 
(Evendale,  Ohio)  under  contract  to  AFAPL/TBC. 

Scope 

Experimental  testing  and  evaluation  of  analytically 
derived  blade  geometries  are  an  integral  part  of  this  AFAPL 
Heat  Transfer  Facility  blade  design  upgrade  program.  How¬ 
ever,  it  is  beyond  the  scope  of  this  thesis  research  effort 
to  conduct  the  experimental  portion  of  the  transonic  blade 
development  program. 

All  flowfield  calculations  and  estimates  are  valid  for 
the  mid-span  2-D  streamsurface  only.  Aerodynamic  performance 
and  flowfield  behavior  are  predicted  for  the  design  point, 
plus  a  number  of  off-design  test  conditions.  The  determin¬ 
ation  of  heat  transfer  loads  and  mechanical  stress  analyses 
are  also  beyond  the  scope  of  this  study. 

The  design  point  environment  can  be  expressed  in  terms 
of  the  following  gas  property  and  velocity  vector  information 
29.0  N/cm2  (42.0  psia)  total  pressure,  13.1  N/cm2  (19.0  psia) 
back  pressure,  1666. 3°K  (3000°R)  total  temperature,  1.26 
specific  heat  ratio,  315.0  m2/sec2-°K  (1883.7  ft2/sec2-°R) 
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gas  constant,  8.0°  gas  incidence  angle,  and  65.0°  gas  exit 
angle. 

Four  geometric  shape  parameters  are  varied  and  opti¬ 
mized  during  the  parametric  design  process:  location  of 
maximum  section  thickness,  stagger  angle,  leading  edge 
bluntness  factor,  and  trailing  edge  taper  factor.  These 
parameters  are  defined  in  Chapter  V.  Performance  sensitivity 
to  off-design  gas  incidence  angle,  total  temperature,  and 
pressure  ratio  are  also  investigated. 

Assumptions 

A  number  of  assumptions  are  inherent  to  the  TDS  flow- 
field  calculation  procedure.  The  most  notable  assumptions 
are:  1)  the  flow  is  steady,  2)  the  flow  is  inviscid,  3)  the 

flow  is  isentropic  (shock-free),  4)  the  fluid  behaves  as 
an  ideal  gas,  and  5)  mass,  momentum  and/or  energy  are  not 
transported  across  a  two-dimensional  streamsurface. 

The  first  assumption  implies  that  the  inlet  velocity 
profile  to  the  test  blade  is  uniform.  Because  the  transition 
section  immediately  upstream  of  the  test  section  is  a  con¬ 
vergent  duct  (see  Figures  1  and  2)  with  favorable  decreasing 
pressure  gradients,  viscous  boundary  layer  growth  has 
negligible  effect  on  the  inlet  velocity  profile.  The  test 
section  is  assumed  to  be  located  far  enough  downstream  that 
any  residual  freestream  turbulance  from  the  upstream  com¬ 
bustion  process  is  damped  out.  A  condition  of  steady,  irro- 
tational  flow  is  assumed  to  exist. 
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The  next  assumption  is  the  shock-free  flowfield.  This 
assumption  breaks  down  near  the  trailing  edge  wake  region  of 
the  cascade.  Viscous  effects  become  significant  in  this 
region  also.  However,  shock-free  isentropic  flow  is  a 
reasonable  assumption  for  the  rest  of  the  cascade  flowfield 
where  velocities  are  subsonic  or  slightly  sonic  and  pressure 
gradients  are  negative. 

Finally,  the  fifth  assumption  is  applied  to  a  theoreti¬ 
cal  blade-to-blade  streamsurf ace  at  the  cascade  mid-span 
location.  This  two  dimensional  flow  assumption  is  justified 
on  the  basis  of  the  rectilinear  test  section  cascade  geometry 
and  the  symmetrically  convergent  upstream  transition  duct. 

The  accuracy  and  validity  of  the  experimental  data 
contained  in  AFAPL-TR-67-147 1  is  tacitly  accepted.  This 
assumption  became  a  necessary  premise  during  the  TDS  veri¬ 
fication  process. 

General  Approach 

The  basic  approach  toward  solving  the  problem  of  de¬ 
signing  a  transonic  turbine  stator  via  analytical  techniques 
is  roughly  a  four  phased  process;  1)  The  flowfield  calcu¬ 
lation  procedure  is  verified  for  correctness;  2)  Blade 
design  criteria  are  developed,  and  the  existing  blade  profile 
deficiencies  are  identified;  3)  Blade  shape  parameters  are 
optimized  along  aerodynamic  considerations  to  yield  a  final, 
transonic  design;  and  4)  Off-design  performance  is  investi¬ 
gated. 


7 


The  first  phase  of  this  design  project  is  concerned 
with  establishing  the  validity  of  the  preceding  initial 
assumptions  and  the  TDS  calculation  techniques.  Verification 
of  the  analytical  calculation  techniques.  Verification  of 
the  analytical  calculation  method  as  well  as  confirmation 
of  the  assumptions  is  achieved  through  correlation  of  sub¬ 
sonic  test  data1  for  the  present  airfoil  cascade  with  TDS 
flowfield  predictions.  Further  substantiation  of  the  TDS 
program  is  obtained  through  comparison  of  TDS  estimates  for 
a  transonic  flow  case  with  General  Electric  Company  test  data, 
and  other  analytical  prediction  programs. 

Accomplishment  of  the  phase  two  objective  requires  the 
development  of  suitable  aerodynamic  performance  criteria. 

Next,  these  criteria  are  applied  to  TDS  generated  flowfield 
predictions  for  the  current  blade  profile  in  a  transonic 
flow  regime  (design  point) . 

Finally,  trends  are  established  in  phase  three  between 
the  effect  of  varying  various  geometric  blade  shape  para¬ 
meters  and  their  corresponding  influence  on  aerodynamic 
performance.  A  final  transonic  design  results  from  opti¬ 
mization  of  the  blade  shape  trends.  Off-design  cases  are 
also  investigated. 

Sequence  of  Presentation 

The  next  chapter,  entitled  "Analytical  Design  Methods", 
contains  a  brief  historical  summary  of  the  development  of 
analytical  calculation  techniques  for  the  turbomachinery 
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cascade  throughflow  problem.  In  addition,  a  description  of 
the  TDS  computer  design  tool  developed  by  General  Electric 
Company  under  U.S.  Air  Force  contract  #F33615-75-C-2073 
appears  in  Chapter  II. 

The  phase  one  results  from  the  TDS  verification  process, 
including  correlation  of  theoretical  predictions  with  ex¬ 
perimental  data  for  the  current  blade  design,  are  presented 
in  Chapter  III. 

Transonic  design  operating  conditions  and  performance 
standards  are  defined  in  the  first  two  sections  of  Chapter  IV. 
These  criteria  are  applied  to  an  analytical  prediction  of  the 
flowfield  surrounding  the  present  design  in  a  transonic  en¬ 
vironment.  The  resulting  deficiencies  and  shortcomings  of 
the  current  blade  are  then  identified  and  documented  in  the 
last  section  of  Chapter  IV. 

The  phase  three  parametric  design  study  results  are 
discussed  in  Chapter  V.  The  final  design  coordinates  appear 
in  both  tabular  and  graphical  form.  Estimates  of  aerodynamic 
performance  sensitivity  to  various  off-design  conditions 
appear  in  Chapter  VI.  The  influence  of  viscous  effects  on 
performance  is  also  addressed. 

Conclusions  and  recommendations  for  further  study  are 
offered  in  Chapter  vn. 
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II.  Analytical  Design  Method 

The  design  of  turbomachine  cascade  sections  has  tradi¬ 
tionally  been  based  upon  empirical  relations  and/or  one¬ 
dimensional  flow  analyses.  Thus,  the  standard  practice  for 
the  design  of  compressor  and  turbine  cascades  has  primarily 
consisted  of  a  "build  and  try"  approach.  The  advent  of 
advanced  gas  turbine  engines  for  aircraft  application,  along 
with  their  staggering  development  costs  and  stringent  perfor¬ 
mance  requirements  (high  temperature  operation,  lightweight, 
high  thermodynamic  efficiency) ,  has  necessitated  a  revolu¬ 
tionary  change  in  compressor  and  turbine  design  procedures. 
Today,  the  performance  of  modern  designs  is  analytically 
predicted  before  a  commitment  is  made  to  proceed  with  an 
expensive  fabrication  and  test  program.  Large,  high  speed, 
digital  computers  have  made  this  approach  more  attractive 
from  a  cost  effectiveness  standpoint. 

Pursuant  to  the  goal  of  analytically  predicting  turbo¬ 
machine  cascade  performance,  a  substantial  number  of 
investigators  from  industry,  government,  and  academia,  have 
developed  cascade  through-flow  calculation  techniques.  There¬ 
fore,  one  objective  of  this  chapter  is  to  present  an  overview 
of  the  various  methods  currently  used  to  compute  turbomachine 
cascade  flowfields.  However,  the  basic  intent  of  this 
section  is  to  provide  a  detailed  description  about  the  nature 
and  operation  of  the  major  analytical  design  tool  (i.e.,  the 
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TDS  computer  programs)  that  will  be  used  to  generate  new 
transonic  blade  profiles  and  to  determine  the  corresponding 
cascade  flowfield  characteristics . 

Cascade  Through-Flow  Calculation  Techniques 

The  cascade  through-flow  problem  is  three-dimensional, 
time-dependent  and  viscous  in  nature.  Most  modern  analysis 
techniques  dealing  with  three-dimensional  cascade  flowfields 
are  based  upon  an  approximate  method  set  forth  by  Wu2  in  1951. 
Wu's  paper  is  primal  to  the  field  because  he  was  one  of  the 
first  investigators  to  make  the  necessary  assumptions  and 
simplifying  approximations  for  transforming  the  overall  three- 
dimensional  problem  into  a  series  of  two  dimensional  compu¬ 
tations  . 

Wu's  basic  through-flow  theory  proposes  that  the  three- 
dimensional  flowfield  should  be  modeled  as  a  set  of  two- 
dimensional  flows;  one  of  the*:e  flows  is  located  in  blade-to- 
blade  surfaces  (the  SI  stream  surfaces) ,  while  the  other 
flow  lies  on  hub-to-shroud  surfaces  (the  S2  stream  surfaces) . 
The  relative  orientation  of  the  Si  and  S2  stream  surfaces  is 
depicted  in  Figure  3. 

Two  classical  solution  approaches  to  the  problem  of 
calculating  flows  on  either  the  SI  or  S2  stream  surfaces  are 
reported  in  literature;  the  matrix  method,  and  the  streamline 
curvature  method.  Iterative  numerical  techniques  are  an  in¬ 
herent  feature  of  both  solution  procedures.  Iterative  methods 
guess  an  initial  stream  function/velocity  distribution  from 
which  the  rest  of  the  fluid  properties  in  the  flowfield  are 
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Figure  3  SI  and  S2  Stream  Surfaces 


determined.  The  equation  of  motion  is  solved  by  either  the 
matrix  or  streamline  curvature  approach,  and  thus,  a  new 
velocity  field  is  predicted.  This  iterative  process  continues 
until  a  set  of  specified  convergence  criteria  are  satisfied. 

Wu2  3  suggested  the  use  of  matrix  methods.  In  this  type 
of  approach,  sometimes  referred  to  as  a  finite-element  method, 
the  iterations  are  performed  on  a  fixed  set  of  grid  lines. 

The  matrix  method  has  recently  been  used  to  compute  cascade 
flowfields  by  Bosman  and  El-Shaarawi4 ,  Bosnian5,  Katsanis  and 
McNally  6  7 ,  and  others8  1 5 .  Gopalakrishnan  and  Bozzola16, 
and  McDonald17  pose  the  through-flow  problem  in  time  dependent 
form  and  use  finite-element  techniques  to  solve  for  transonic 
flows  in  turbomachinery  cascades.  Ives  and  Liutermoza18  use 
conformal  mapping  techniques  in  conjunction  with  a  finite- 
difference  grid  to  obtain  transonic  cascade  solutions. 

In  the  streamline  curvature  method,  the  grid  is  composed 
of  the  streamlines  themselves  and  a  set  of  semi-orthogonals . 
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During  the  iteration  sequence,  the  grid  lines  and  points  are 
repositioned  after  each  iteration  to  satisfy  mass,  momentum 
and  energy  balances.  Novak  and  Hearsey19,  Smith20,  Frost21, 
and  Wilkinson22  have  reported  the  development  of  computer 
programs  which  employ  streamline  curvature  techniques. 

The  concept  of  quasi-three-dimensional,  numerical 
solutions  of  flow  through  cascades  refers  to  the  simultaneous 
solution  of  flow  on  the  SI  and  S2  streamsheets .  A  number 
of  investigators4,19  use  coordinated  computer  programs  in 
which  the  SI  stream  surface  shape  is  dependent  on  the  S2 
solutions  while  the  S2  shape  is  dependent  on  the  SI  solutions. 
These  methods  are  indicative  of  the  state-of-the-art  and  have 
established  good  convergence  for  steady,  inviscid  and  sub¬ 
sonic  flow  conditions. 

Turbine  Design  System  (TPS) 

The  specific  intent  of  this  research  is  to  establish  an 
analytical  procedure  for  designing  transonic  turbine  blades. 
Fundamental  geometric  blade  shape  parameters  were  identified, 
and  the  influence  of  these  parameters  upon  cascade  performance 
for  a  given  set  of  operating  conditions  was  determined.  The 
final  design  is  simply  the  geometric  parameter  set  which 
yields  the  highest  amount  of  aerodynamic  performance. 

A  parametric  variation  and  evaluation  study  was  conducted 
to  realize  the  basic  intent  outlined  above.  The  major 
analytical  tool  that  was  used  to  accomplish  this  study  is  the 
Turbine  Design  System  (TDS)  computer  programs  developed  by 
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General  Electric  Company  under  contract  to  AFAPL/Turbine 
Components  Branch.  Two  TDS  subprograms  were  utilized: 
the  BLDDES  and  CASC  modules . 

Given  a  relatively  simple  blade  shape  definition,  the 
Blade  Design  (BLDDES)  module  generates  an  output  file  with 
a  standardized  set  of  blade  shape  coordinates  and  surface 
curvature  information.  The  Cascade  Analysis  by  Streamline 
Curvature  (CASC)  module  returns  the  inviscid,  shock-free 
flowfield  solution  for  a  cascade  of  airfoils  defined  by  the 
BLDDES  output  file.  Figure  4  is  a  schematic  diagram  which 
graphically  displays  the  function  and  role  of  each  module. 
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Further  information  about  the  operation  of  the  TDS  computer 
programs  can  be  found  in  the  TDS  User's  Manual23  and  AFAPL 
TR-78-9224. 

BLDDES  MODULE.  Two  modes  of  data  input  tc  the  BLDDES 
module  are  possible.  Input  may  consist  of  a  parametric 
blade  definition  or  as  an  array  of  digitized  blade  coordinate 
data  taken  from  a  layout  drawing  or  design  report.  The  latter 
mode  of  data  input  was  used  to  obtain  a  standardized  defini¬ 
tion  of  the  current  test  blade  profile.  This  information  was 
required  to  perform  the  TDS  verification  process  (Chapter  III) 
and  to  arrive  at  an  analytic  evaluation  of  the  current  blade 
at  transonic  operating  conditions  (Chapter  IV) .  The  para¬ 
metric  definition  input  mode  option  was  utilized  to  generate 
blade  profiles  for  aerodynamic  analysis  during  the  parametric 
design  study  (Chapter  V) . 

The  standard  definition  returned  by  BLDDES  is  composed 
of  blade  surface  coordinates  that  are  tabulated  with  respect 
to  the  blade  section  stacking  point.  BLDDES  output  is  con¬ 
sistent  with  the  input  requirements  for  the  aerodynamic 
analysis  module;  CASC. 

CASC  Module.  The  CASC  module  yields  steady-state,  in- 
viscid,  shock-free  flow  solutions  for  a  cascade  of  airfoils. 
CASC  solutions  for  any  particular  blade  section  will  lie  on 
the  blade-to-blade  SI  streamsurf ace  normal  to  the  blade 
(see  Figure  3) . 

The  streamline  curvature  calculation  technique  employed 
by  CASC  requires  the  input  of  accurately  determined  boundary 
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surface  coordinates  from  the  standardized  BLDDES  output  $BB 
Namelist.  The  position  and  curvature  of  the  streamlines  are 
found  such  that  the  continuity  equation,  the  cross-stream 
momentum  equations,  and  the  boundary  surface  curvature  in¬ 
formation  are  satisfied. 

Although  the  inviscid  assumption  costs  CASC  a  small 
amount  of  inaccuracy,  it  reduces  the  complexity  of  the  cal¬ 
culations  considerably.  CASC's  main  asset  is  that  it  requires 
substantially  less  central  processor  time  and  core  space  to 
process  a  given  airfoil  contour  as  compared  with  other  aero¬ 
dynamic  analysis  programs  reported  in  literature.  Thus,  the 
"quick-look"  capability  offered  by  CASC  makes  the  program 
ideally  suited  towards  performing  a  parametric  design  study 
involving  a  number  of  cases  to  be  analyzed. 
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III.  TPS  Verification  Investigation 


Purpose 

Verification  of  the  General  Electric  TDS  computer  pro¬ 
grams  was  included  as  a  major  step  in  this  research  study 
for  two  reasons.  First,  the  TDS  codes  had  been  recently 
procured  by  the  U.S.  Air  Force  in  November  1977,  and  had  not 
been  sufficiently  verified  for  accuracy  and  consistency. 
Second,  it  was  necessary  to  determine  whether  or  not  the 
Heat  Transfer  Facility's  rectalinear  cascade  geometry  along 
with  the  simplifying  inlet  velocity  profile  assumptions  could 
be  analyzed  by  TDS  to  yield  accurate  results. 

Pressure  distributions  for  the  existing  turbine  airfoil 
were  experimentally  determined  by  Quick,  Henderson,  and  Tall 
during  a  1967  cold-flow  investigation  study.  The  results 
are  tabulated  in  AFAPL  Technical  Report  67-147 l. 

Hence,  the  overall  TDS  verification  objective  was  to 
obtain  pressure  distribution  predictions  for  the  present 
blade  geometry  that  match  available  experimental  data  at  two 
different  sets  of  operating  conditions.  Comparison  of  analy¬ 
tical  TDS  pressure  distribution  predictions  with  actual  Heat 
Transfer  Facility  experimental  data  was  expected  to  substan¬ 
tiate  the  validity  of  the  initial  assumptions  and  the  computer 
programs . 
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Procedure 


The  general  approach  used  to  verify  the  TDS  programs 
consisted  of  three  distinct  phases.  First,  thermodynamic 
and  geometric  boundary  conditions,  as  well  as  an  inlet 
velocity  profile,  were  specified  for  input  into  the  TDS 
programs.  These  boundary  conditions  were  chosen  so  that 
they  were  a  direct  reflection  of  the  test  conditions  under 
which  the  experimental  data  were  obtained. 

Next,  verification  criteria  were  developed.  In  theory, 
the  analytical  codes  should  yield  matching  pressure  distri¬ 
butions  for  the  two  known  sets  of  test  conditions  that  were 
considered.  Additionally,  the  computer  programs  should 
yield  predictable  flowfiald  solutions  for  small  variations  of 
the  inlet  boundary  conditions. 

Finally,  the  results  of  the  verification  investigation 
were  evaluated.  The  conclusions  are  presented  in  the  last 
section  of  this  chapter. 

Boundary  Conditions 

Thermodynamic .  During  the  1967  experimental  cold  flow 
investigation  of  the  present  turbine  airfoil1,  pressure 
distributions  were  obtained  for  various  gas  operating  con¬ 
ditions.  Two  of  these  cases  were  selected  for  analytical 
correlation  purposes.  The  thermodynamic  properties  for  each 
of  the  two  cases  are  tabluated  in  Tables  I  and  II. 

Geometric.  Establishment  of  geometric  boundary  condi¬ 
tions  for  aerodynamic  analysis  by  the  CASC  module  consists 


of  defining  the  current  blade  profile,  and  then  standardizing 
the  definition  with  respect  to  a  stacking  point  origin  via 
the  BLDDES  module.  The  portion  of  the  BLDDES  input  $BB 
namelist  (FORTRAN)  that  defines  the  present  turbine  airfoil 
shape  appears  in  Table  III.  Figure  5  is  a  scaled  drawing 
of  the  existing  turbine  airfoil  section. 

The  philosophy  used  to  define  the  current  blade  section 
is  the  following;  1)  represent  the  leading  edge  as  an  ellipse, 
2)  represent  the  trailing  edge  as  a  circle,  and  3)  input  the 
pressure  and  suction  surfaces  as  digitized  coordinate  pairs. 
This  philosophy  becomes  obvious  when  the  quantities  contained 
in  Table  III  are  clarified  in  the  following  paragraphs. 

The  pressure  and  suction  surfaces  are  defined  by  the 
(X, Y)  pairs  in  the  PRESS  and  SUCTN  arrays.  These  coordinates 
are  taken  from  the  original  Allison  Division  blade  design 
report25.  The  coordinate  system  that  was  employed  fixes  the 
X-axis  parallel  to  the  axial  centerline,  positions  the  leading 
edge  highlight  point  at  (0,3),  and  the  suction  surface  is  the 
upper  surface . 

ISTK  is  the  stacking  point  indicator.  A  value  equal 
to  unity  locates  the  stacking  point  at  the  center  of  gravity 
of  the  airfoil  section.  The  choice  of  stacking  point  location 
is  academic;  since  choice  of  any  stacking  point  location  will 
yield  identical  stacking  results  for  a  rectalinear  two- 
dimensional  cascade. 

A  zero  value  for  IVANE  notifies  the  software  logic  that 
the  airfoil  has  a  blade  input  orientation,  i.e.,  the  suction 
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Table  III 


Geometric  Boundary  Conditions 


PRESS= 


0.0740 

2.9205 

0.1273 

2.9126 

0.1826 

0.2337 

2.8945 

0.2851 

2.8822 

0.3342 

0.3827 

2.8558 

0.4305 

2.8404 

0.4771 

0.5232 

2.8045 

0.5912 

2.7729 

0.6574 

0.7399 

2.6859 

0.8153 

2.6334 

0.8910 

0.9654 

2.5131 

1.0348 

2.4479 

1.1023 

1.1673 

2.3093 

1.2301 

2.2357 

1.2879 

1.3429 

2.0802 

1.4015 

1.9865 

1.4317 

1.4619 

1.8815 

SUCTN= 

0.0466 

3.0730 

0.1016 

3.0984 

0.1665 

0.2294 

3.1380 

0.3027 

3.1420 

0.3728 

0.4422 

3.1343 

0.5108 

3.1210 

0.5787 

0.6442 

3.0741 

0.7366 

3.0228 

0.8203 

0.9053 

2.8744 

0.9709 

2.7994 

1.0362 

1.1006 

2.6381 

1.1611 

2.5568 

1.2194 

1.2746 

2.3894 

1.3274 

2.3034 

1.3771 

1.4253 

2.1294 

1.4781 

2.0278 

1.5054 

1.5327 

1.9209 

ISTK=1 
IVANE=0 
LEDGE=1 
ITMODE=l 
TE=0. 08335 
E=1 . 45 
EPSI=8 . 0 
SCALIN=1.0 
ZLE=0 . 0 
ZTE=1 . 5327 
RLE=23. 87324 
RREF=23. 87324 
RTE=*23. 87324 
NBLDS=100 


2.9033 

2.8698 

2.8234 

2.7372 

2.5756 

2.3801 

2.1593 

1.9340 


3.1205 

3.1418 

3.1006 

2.9590 

2.7195 

2.4742 

2.2167 

1.9744 
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Figure  5  Existing  Turbine  Airfoil  Section 


side  is  uppermost.  LEDGE  equal  to  one  means  that  the  PRESS 
and  SUCTN  coordinate  arrays  do  not  include  the  leading  and 
trailing  edge.  BLDDES  will  then  generate  an  ellipse  for  the 
leading  edge  while  either  a  square  or  circle  will  be  fitted 
to  the  trailing  edge. 

ITMODE  is  the  trailing  edge  shape  indicator.  A  value 
of  one  causes  the  BLDDES  programming  to  define  the  trailing 
edge  as  a  circle  with  diameter  equal  to  the  trailing  edge 
thickness  (TE) .  From  a  plot  of  the  original  design  coordi¬ 
nates,  the  trailing  edge  thickness,  TE,  was  determined  to 
be  0.2117  cm  (0.08335  in). 

The  eccentricity  parameter,  E,  is  the  ratio  of  leading 
edge  ellipse  semi-major  to  semi -minor  axes.  After  a  plot  of 
the  blade  was  made  from  coordinate  data  and  several  iter- 
tions  were  performed  to  obtain  the  best  curve  fit,  a  value 
for  E  of  1.45  was  determined.  EPSI  is  an  estimate  of  the 
leading  edge  ellipse  semi -major  axis  inclination  angle. 

SCALIN  defines  the  scale  size  of  the  input  data.  Actual 
coordinates  were  used.  ZLE  and  ZTE  refer  to  the  axial  location 
of  the  leading  and  trailing  edge  highlight  planes.  RLE, 

RREF,  and  RTE  are  the  radii  which  define  the  airfoil  section 
radial  location  in  a  hypothetical  curvalinear  engine.  Because 
the  facility  test  section  is  rectalinear,  a  large  radius 
value  was  chosen  to  simulate  this  condition.  Thus,  a  radius 
value  of  60.638  cm  (23.87324  in)  was  chosen  along  with  100 
blades  (NBLDS)  to  yield  a  pitch  value  of  3.81  cm  (1.5  in) 
which  is  the  blade  spacing  in  the  two-dimensional  cascade. 
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Inlet  Velocity  Profile .  A  condition  of  steady,  turbu¬ 
lent,  and  fully-developed  flow  has  been  assumed  to  exist  at 
the  test  section  inlet.  Because  fully-developed,  turbulent 
velocity  profiles  are  characterized  by  fairly  uniform  core 
regions,  the  test  blade  inlet  velocity  profile  is  a  set  of 
constant,  uniform  velocity  vectors.  The  direction  of  the 
velocity  vectors  is  the  gas  incidence  angle,  Py  ,  and  the 
magnitude  is  expressed  in  non-dimensional  form  as  the  inlet 
Mach  number,  Mx . 

Due  to  the  geometry  of  the  upstream  transition  section, 
a  slight  gas  incidence  angle  is  known  to  exist.  The  gas 
incidence  angle  is  estimated  to  be  8°  +  8°.  As  will  be 
shown  in  the  theoretical  development  of  the  third  verifica¬ 
tion  criterion,  small  gas  incidence  angles  (+15°)  have 
negligible  effect  on  the  pressure  distribution  of  an  airfoil 
in  a  turbine  cascade.  Therefore,  the  amount  of  error  as¬ 
sociated  with  +8°  is  insignificant.  Because  gas  incidence 
angle  is  primarily  a  function  of  hardware  geometry  and  not 
gas  operating  conditions,  it  is  assumed  to  remain  constant 
for  both  experimental  cases  under  consideration. 

The  inlet  Mach  number  was  approximated  for  each  case 
through  application  of  the  following  one-dimensional,  isen- 
tropic  equation; 
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where  ift  is  the  mass  flow  rate,  A  is  the  cross-sectional 
area  normal  to  the  flow,  T0  is  the  total  temperature,  R  is 
the  universal  gas  constant  divided  by  the  molecular  weight 
of  air,  y  is  the  specific  heat  ratio,  and  PQ  is  the  total 
pressure.  When  the  appropriate  Table  I  and  II  values  are 
inserted  into  equation  one,  along  with  a  cross-sectional  area 
of  148.4  cm2  (23  in2),  the  inlet  Mach  number  is  estimated 
to  be  0.29  and  0.30  for  cases  one  and  two  respectively. 

Because  of  the  significant  amount  of  error  associated  with 
the  figures  for  mass  flow  rate  and  the  cross-sectional  area 
normal  to  the  flow,  these  approximations  are  only  useful  as 
a  general  estimate  of  the  inlet  Mach  number  level.  Therefore, 
a  number  of  cases  with  varying  inlet  mach  numbers  were  anal¬ 
yzed  before  a  precise  analytical-experimental  correlation  was 
reached. 

The  case  one  and  two  BLDDES  output  files  (CASC  input 
files)  which  contain  the  standardized  blade  definition  of 
the  existing  airfoil,  plus  the  respective  thermodynamic 
and  inlet  velocity  profile  boundary  information,  are  listed 
in  Appendix  A.  Further  explanation  of  the  Appendix  A  $BB 
namelist  variables  can  be  found  in  the  TDS  User's  Manual23. 

Standards 

Criterion  One .  For  a  given  input  combination  of  thermo¬ 
dynamic,  geometric  and  inlet  velocity  profile  information, 
the  CASC  module  returns  a  value  for  the  cascade  exit  static 
pressure,  P2 ,  the  gas  exit  angle,  )S2,  plus  a  listing  of  the 
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airfoil  pressure  distribution.  For  each  of  the  two  test 
cases,  the  first  verification  criterion  requires  that  the 
analytical  pressure  distribution  prediction  matches  the  ex¬ 
perimental  data.  Analytical-experimental  correlation  must 
occur  for  the  boundary  conditions  stated  in  the  previous 
section,  plus  a  small  change  in  the  inlet  Mach  number  value. 

Since  the  exit  static  pressure,  P2,  must  equal  the 
cascade  back  pressure,  P^,  for  a  subsonic  cascade  operating 
condition,  the  exit  static  pressure  is  equal  to  the  atmos¬ 
pheric  pressure,  plus  exhaust  duct  line  losses,  ^P; 

P2  =  Pb  =  Pa  +  (2) 

Because  case  one  is  at  a  lower  mass  flow  rate  operating 
condition  than  is  case  two,  lower  exhaust  duct  exit  veloci¬ 
ties  will  exist  during  case  one.  Thus,  the  line  pressure 
losses  are  expected  to  be  less  for  case  one.  In  addition, 
the  line  pressure  losses  for  either  case  should  be  greater 
than  zero,  but  less  than  20  percent  of  the  atmospheric 
pressure  level.  Therefore,  another  requirement  of  the  first 
verification  criterion  is  that  the  following  inequality  must 
be  staisfied; 


0.0  <  4PCase  one  ^ 4pcase  two  ^  0*20  Pa  (3) 

Criterion  Two .  The  gas  exit  angle,  /J2 »  is  primarily 
fixed  by  the  cascade  blade  outlet  angle,  Cf2.  The  blade 
outlet  angle  is  defined  in  Figure  6  as  the  angle  formed  by 
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the  tangent  line  to  the  trailing  edge  camberline  and  the 
axial  chord  line.  The  existing  turbine  cascade  has  an  out¬ 
let  blade  angle  of  61.2°. 

In  most  turbine  and  compressor  cascades,  a  slight  gas 
deviation  angle,  d,  will  exist.  Gas  deviation  angles  for 
subsonic  cascades  are  small  and  result  from  the  pressure 
equilibrium  process  that  occurs  in  the  trailing  edge  wake 
region  of  the  blades.  Vincent26  estimates  that  deviation 
angles  are  +2°  for  most  turbine  cascades.  Cox27  reports 
deviation  angles  of  the  order  of  5°  from  cascade  test  data. 

Hence,  the  second  verification  criterion  requires  that 
the  gas  exit  angle  computed  by  CASC  be  equal  to  61.2°  +  6, 
where : 

-5°  <  6  <  5°  (4) 


Criterion  Three.  The  governing  two-dimensional  conser¬ 
vation  of  mass  law  is  expressed  in  differential  form  by  the 
following  equation: 


dp  +  dpu  dpv 
3*t  dx  dy 


0 


(5) 


When  the  preceeding  continuity  equation  is  applied  to 
the  mid-span  turbine  cascade  control  volume  depicted  in 
Figure  6,  the  equation  reduces  to  the  following  steady-state, 
one-dimensional  expression: 
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Equation  7  is  a  simple  continuity  balance  which  states 
that  the  mass  flow  through  the  inlet  plane  is  equal  to  mass 
flow  out  the  exit  plane.  Upon  substitution  of  isentropic 
gas  relationships  and  an  algerbraic  rearrangement  of  terms 
(see  Appendix  B) ,  Equation  7  becomes; 
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where 

f3 1  =  gas  incidence  angle 
Mi  =  inlet  Mach  number 
/?2  =  gas  exit  angle 
M2  =  exit  Mach  number 

Equation  8  is  non-dimensional  and  it  relates  the  cascade 
inlet  conditions  to  the  exit  conditions. 

The  pressure  distribution  around  an  airfoil  of  a  typical 
turbine  cascade  has  been  observed  to  be  a  function  of  chan¬ 
nel  flow  characteristics  rather  than  parameters  associated 
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with  external  airfoil  theory.  Moreover,  for  a  given  airfoil 
in  a  cascade,  the  pressure  coefficient  distribution  is  a 
function  of  one  unique  cascade  pressure  ratio  (i.e.,  mass 
flow  rate) : 


Cp  Distribution  =  f (PR) 


(9) 


where 


Cp  =  pressure  coefficient 
PR  =  cascade  pressure  ratio 


The  pressure  coefficient,  Cp,  is  defined  as  the  local 
static  pressure  divided  by  the  total  pressure.  Cascade 
pressure  ratio,  PR,  is  equal  to  the  exit  static  pressure 
divided  by  the  total  pressure.  Exit  Mach  number,  M2,  can  be 
defined  in  terms  of  the  cascade  pressure  ratio,  P2/P0/  by 
the  following  relationship: 


or 


(10) 


(11) 


where 


y  =  1.4 
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Two  points  regarding  Equation  8  require  further  clari¬ 
fication.  First,  the  right  hand  side  of  Equation  8  remains 
constant  for  a  given  pressure  distribution  (matching  or 
otherwise)  because  the  gas  exit  angle  is  relatively  invari¬ 
ant  with  varying  inlet  combinations  and  the  exit  Mach  number 
is  also  constant  as  defined  by  Equations  11  and  9.  Second, 
for  small  gas  incidence  angles ,  cos  ,  is  approximately 
equal  to  unity.  Thus,  turbine  cascades  are  insensitive  to 
small  gas  incidence  angles  and  the  inlet  Mach  number  is 
relatively  fixed  for  a  given  pressure  distribution. 

Consequently,  if  a  matching  pressure  distribution  is 
found  such  that  the  first  two  verification  criteria  are 
satisfied,  the  gas  exit  angle,  02r  and  the  cascade  pressure 
ratio,  P2/P0,  returned  by  CASC  are  a  pair  of  unique  quantities. 
Substitution  of  these  gas  exit  conditions  into  Equations  11 
and  then  8  will  yield  a  relationship  which  involves  only 
two  variables;  the  inlet  Mach  number  and  the  gas  incidence 
angle.  Hence,  the  third  verification  criterion  requires 
that  the  TDS  programs  should  yield  matching  pressure  distri¬ 
butions  for  the  inlet  velocity  profile  combinations  defined 
by  Equation  8,  where  the  right  side  of  the  equation  is  a 
previously  determined  constant. 

Results 

For  the  specified  boundary  conditions  which  define  either 
verification  case,  analytical  pressure  distributions  were 
obtained  that  match  the  experimental  data.  At  a  gas  incidence 
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angle  of  eight  degrees ,  the  inlet  Mach  number  values  which 
yield  the  matching  data  are  0.2842  and  0.2990  respectively. 
Note  that  the  difference  between  these  inlet  Mach  number  val¬ 
ues  and  the  Equation  1  approximations  are  negligible. 

The  CASC  output  files  for  both  cases  are  presented  in 
Appendix  C.  Figures  7  and  8  are  pressure  coefficient  dis¬ 
tribution  plots  which  graphically  demonstrate  the  correlation 
between  the  experimental  and  analytical  results. 

Comparison  of  the  Figure  7  and  8  curves  indicates  that  a 
small  (1-5  percent)  discrepancy  exists  between  predicted  and 
and  actual  data.  The  analytic  pressure  coefficients  are 
slightly  higher  than  the  corresponding  experimental  results. 
This  situation  occurs  primarily  because  CASC  is  an  inviscid 
analysis  tool.  Due  to  the  fact  that  CASC  neglects  the  pres¬ 
ence  of  boundary  layers,  it  does  not  consider  displacement 
thickness  effects.  Thus,  lower  velocities  (higher  pressures) 
are  predicted  due  to  larger  effective  flow  areas. 

The  six  to  eight  percent  error  associated  with  the  pre¬ 
diction  of  the  minimum  throat  location  is  a  function  of  three 
probable  causes.  First,  boundary  layer  induced  displacement 
thickness  can  cause  the  minimum  throat  location  to  shift  for¬ 
ward.  Second,  local  blade  surface  curvature  effects  from 
manufacture  of  the  test  specimen  could  also  have  advanced  the 
minimum  area  location.  Finally,  small  inaccuracies  inherent 
in  the  experimental  procedure  could  have  accounted  for  mis¬ 
representation  of  the  minimum  throat  location. 
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Figure  7  Case  One  Pressure  Coefficient  Distribution 
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Figure  8  Case  Two  Pressure  Coefficient  Distribution 
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Table  IV 


Trailing  Edge  Gas  Properties 


Test 

Case 

P2 

(N/cm2 ) 

(psia) 

m2 

02 

(deg) 

1 

10.434 

15.1 

0.7053 

59.64 

2 

11.052 

16.0 

0.8275 

60.08 

Table  IV  contains  the  trailing  edge  conditions  returned 
by  CASC  for  the  matching  pressure  coefficient  distributions 
plotted  in  Figures  7  and  8.  Line  pressure  losses,  Ap,  are 
calculated  by  inserting  the  Table  IV  exit  static  pressures 
into  Equation  2.  For  an  atmospheric  pressure  of  9.885  N/ 
cm2  (14.3  psia) ,  the  line  losses  are  0.549  N/cm2  (0.8  psi) 
and  1.167  N/cm2  (1.7  psi).  These  4p  values  satisfy  the 
relationship  expressed  by  Equation  3. 

The  gas  deviation  angle,  <5,  is  the  difference  between 
the  blade  outlet  angle  and  the  gas  exit  angle.  The  computed 
deviation  angles  of  1.56°  and  1.12°  for  the  two  test  cases 
are  within  the  Equation  4  requirements. 

In  order  to  satisfy  the  third  verification  criterion, 
the  effect  of  varying  the  inlet  Mach  number  and  gas  incidence 
angle  was  investigated.  A  combined  total  of  ninety-one  dif¬ 
ferent  inlet  combinations  were  analyzed  for  both  cases. 
Appendicies  D  and  E  contain  the  pressure  coefficient  distri- 
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Figure  9  Case  One  Inlet  Velocity  Profiles 

Yielding  Matching  Cp  Distributions 
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bution  plots  for  a  representative  sampling  of  the  ninety- 
one  point  data  field. 

Figures  9  and  10  were  generated  from  a  visual  inspection 
of  the  entire  set  of  pressure  coefficient  distributions. 
Distributions  which  match  the  appropriate  set  of  experimental 
data  are  plotted  as  points  on  a  versus  Mi  coordinate  system. 
The  matching  data  collapse  along  an  exit  static  pressure  (P2) 
curve  whose  value  is  listed  in  Table  IV.  This  result  is  pre¬ 
dicted  by  the  Equation  9  functional  relationship. 

A  line  defined  by  Equation  8,  where  the  right  side  of  the 
equation  has  been  evaluated  by  the  Table  IV  and  M2  con¬ 
straints,  is  superimposed  upon  the  matching  data  points. 

The  agreement  between  the  locus  of  expected  and  calculated 
matching  pressure  distributions  is  excellent.  Note  that  the 
influence  of  small  gas  incidence  angles  upon  the  pressure 
distribution  is  negligible,  while  slight  Mach  number  changes 
will  define  totally  different  pressure  distribution  curves. 

Conclusions 

The  results  presented  in  the  previous  section  satisfy 
all  three  verification  criteria  for  the  two  test  cases  con¬ 
sidered.  In  addition,  the  simplifying  inlet  velocity  profile 
assumptions  in  the  test  facility  have  been  substantiated. 

Thus,  the  ability  of  the  TDS  computer  programs  to  calculate 
two-dimensional  cascade  flowfields  accurately  in  the  Heat 
Transfer  Facility  test  environment  has  been  established. 
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IV.  Blade  Design  Criteria 

The  objective  of  this  blade  design  study  is  to  provide 
the  Air  Force  Aero  Propulsion  Laboratory  with  a  updated  blade 
profile  to  be  used  in  an  "in-house"  turbine  cooling  test 
program.  The  final  blade  design  must  be  compatible  with  the 
static  test  facility  in  which  the  blades  are  to  be  tested. 

The  airfoil  shape  must  also  exhibit  maximum  aerodynamic  per¬ 
formance  in  the  transonic  flow  regime  specified  by  AFAPL/TBC. 

In  order  to  accomplish  the  objectives  stated  above, 
certain  criteria  must  be  established  and  incorporated  into 
the  design  process.  The  threefold  intent  of  this  chapter  is 
to  define  the  specified  transonic  design  point,  to  present 
the  aerodynamic  performance  standards  used  for  selecting 
superior  designs,  and  to  document  a  transonic  evaluation  of 
the  existing  design  using  the  previously  identified  perfor¬ 
mance  standards. 

Design  Point  Operating  Conditions 

The  primary  conditions  which  governed  the  design  of  the 
turbine  airfoil  were  based  on  the  anticipated  operating  con¬ 
ditions  of  the  test  facility  in  which  the  blades  are  to  be 
tested.  A  compilation  of  the  design  point  gas  properties 
specified  by  AFAPL/TBC  appears  in  Table  V. 

The  total  pressure  and  back  pressure  values  listed  in 
Table  V  are  within  the  operating  limits  of  the  test  facility 
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for  a  mass  flow  rate  less  than  2.27  kgm/sec  (5  lbm/sec) . 
Together,  the  two  quantities  define  a  transonic  environment 
with  an  exit  Mach  number  equal  to  1.169. 

The  design  point  total  gas  temperature  of  1666. 3°K 
(3000°R)  is  substantially  greater  than  most  modern  turbine 
inlet  temperatures.  However,  advanced  gas  turbines  are  ex¬ 
pected  to  encounter  this  type  of  temperature  environment. 

The  AFAPL  Heat  Transfer  Facility  has  a  hot  gas  temperature 
generation  capability  of  2500°K  (4500°R)  for  stoichometric 
hydrogen-oxidant  ratios. 

The  mass  flow  rate  limitation  of  2.27  kgm/sec  (5  lbm/sec) 
is  imposed  as  a  constraint  because  this  is  the  maximum  mass 
flow  rate  that  the  facility's  compressors  can  sustain  while 
maintaining  a  total  pressure  level  greater  than  29.0  N/ 
cm2  (42.0  psia) .  The  gas  constant  and  specific  heat  ratio 
are  derived  from  the  hydrogen  combustion  reaction  which  pro¬ 
duces  a  cascade  total  inlet  temperature  of  1666. 3°K  (3000°R) . 

Although  small  gas  incidence  angles  have  negligible  ef¬ 
fect  on  turbine  cascade  performance  (determined  in  Chapter 
III) ,  a  design  point  gas  incidence  angle  of  eight  degrees  was 
chosen  as  a  best  estimate  in  lieu  of  firm  evidence  to  the 
contrary.  When  the  maximum  mass  flow  rate,  along  with  the 
other  Table  V  thermodynamic  quantities,  are  substituted  into 
Equation  1,  the  maximum  allowable  design  point  inlet  Mach 
number  is  found  to  be  0.46. 


Aerodynamic  Performance  Standards 

Transonic  blade  performance  is  usually  measured  exper¬ 
imentally  because  of  the  inability  of  present  day  analytical 
procedures  to  successfully  model  complicated  blade  trailing 
edge  flowfields  where  most  pressure  losses  occur.  Streamline 
curvature  flow  calculation  programs,  including  CASC ,  require 
simplifying  shock-free  and  inviscid  assumptions  which  make 
it  impossible  to  compute  an  estimate  of  total  pressure  losses. 
Thus,  the  establishment  of  aerodynamic  performance  standards 
based  on  minimum  pressure  loss  gives  rise  to  a  fundamental 
question;  can  the  standards  be  applied  to  isentropic  loss- 
free  CASC  predictions?  The  purpose  of  this  section  is  to 
answer  the  question  by  presenting  a  generalized  discussion 
of  the  theory  behind  the  evaluation  of  transonic  criteria 
that  will  be  applied  to  select  superior  designs. 

Trailing  Edge  Pressure  Gradient .  The  primary  function 
of  turbine  stater  blade  rows  is  to  accelerate  the  working 
gas  medium  across  a  static  pressure  drop.  For  a  turbine 
cascade  where  the  flow  is  subsonic  (low  pressure  ratios) ,  the 
fluid  is  accelerated  through  a  convergent  flow  passage  as 
shown  in  Figures  2  and  6.  The  total  pressure  loss  in  a  sub¬ 
sonic  turbine  stator  row  is  very  small.  Viscous  losses  are 
of  the  order  of  five  percent.  There  are  no  shock-related 
losses  because  the  flow  is  everywhere  subsonic.  Finally, 
losses  due  to  separation  usually  do  not  occur  because  the 
pressure  gradients  in  an  accelerating  flowfield  are  favorable. 
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If  the  pressure  ratio  across  a  blade  row  is  large  enough 
to  choke  the  flow,  a  sonic  line  will  appear  at  the  flow  pas¬ 
sage  throat  location.  In  a  typical  convergent  turbine  cascade, 
the  throat  (sonic  line)  is  formed  at  the  trailing  edge  of  a 
blade’s  pressure  surface  and  extends  to  the  suction  surface 
of  the  adjacent  blade.  The  transonic  flowfield  surrounding 
the  trailing  edge  of  a  blade  in  a  convergent  turbine  stator 
row  is  depicted  in  Figure  11. 

At  the  trailing  edge  of  the  Figure  11  blades,  the  flow 
over  the  suction  surface  is  supersonic  (M>1).  The  flow  over 
the  pressure  surface  is  subsonic  (M<1)  until  a  sonic  condi¬ 
tion  is  realized  at  the  trailing  edge.  Thus,  a  region  of 
relatively  high  pressure  exists  on  the  lower  surface,  while 
a  region  of  low  pressure  is  present  on  the  upper  surface. 

Hence,  the  existence  of  an  extremely  large  pressure  gradient 
at  the  trailing  edge  of  the  blades  would  seem  to  be  indicated. 

According  to  Amana,  Demuren  and  Louis28,  this  condition 
is  reconciled  by  an  expansion  of  the  flow  around  the  pressure 
surface  trailing  edge  and  then  a  subsequent  compression  as 
the  flow  turns  into  the  downstream  wake  (see  Figure  11) .  This 
phenomena  occurs  in  order  to  satisfy  the  Kutta  condition  of 
zero  pressure  gradient  at  the  trailing  edge.  The  compression 
shock  wave  formed  during  this  process  originates  from  the 
trailing  edge  wake  region  and  then  impinges  on  the  suction 
surface  of  the  adjacent  blade.  The  boundary  layer-shock 
interation  on  the  suction  surface  of  the  blade  usually 


Figure  11  Trailing  Edge  Transonic  Flowfield 


results  in  flow  separation  or  a  change  in  the  boundary  layer 
thickness.  This,  in  turn,  changes  the  suction  surface  trail¬ 
ing  edge  boundary  conditions  and  thus,  fluctuating  (time- 
dependent)  trailing  edge  pressure  gradients  are  set  up.  This 
is  the  loss  mechanism  which  substantially  reduces  the  perfor-  ' 
mance  of  convergent  cascades  under  transonic  pressure  ratios. 
Amana  et  al28  and  Heinemann29  have  reported  this  phenomenon 
in  literature  and  they  also  present  corroborating  photographic 
evidence. 

Although  the  CASC  program  is  unable  to  predict  the  total 
pressure  losses  which  accompany  these  complicated  transonic 
flow  patterns,  inviscid,  shock-free  theory  can  yield  a  rough 
estimate  of  the  pressure  gradient  that  must  be  reconciled 
at  the  blade  trailing  edge  due  to  poor  design.  Because  the 
pressure  gradient  is  primarily  caused  by  isentropic  flow  over¬ 
expansion  on  the  suction  surface  and  under expans ion  on  the 
pressure  surface,  the  CASC  estimate  of  pressure  gradient  mag¬ 
nitude  is  an  indication  of  the  magnitude  of  the  resultant 
flow  turning,  shock  strength,  and  total  pressure  losses  that 
must  occur  to  satisfy  the  Kutta  condition.  Therefore,  the 
first  aerodynamic  performance  criterion  requires  that  the 
predicted  trailing  edge  pressure  gradient  is  zero. 

i)  -  °  <12> 

/  T.E. 


In  order  to  satisfy  the  criterion  represented  by  equa¬ 
tion  12,  the  flow  on  either  surface  must  be  expanded  to  a 
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single  static  pressure,  i.e.,  the  blade  row  back  pressure. 
Thus,  the  flow  passage  through  the  stator  row  will  necessar¬ 
ily  be  convergent-divergent  in  nature  for  transonic  pressure 
ratios . 

Suction  Surface  Diffusion.  Any  number  of  blade  contours 
can  be  designed  to  satisfy  the  first  performance  criterion 
at  the  design  point  described  in  the  first  section  of  this 
report.  Thus,  a  second  performance  criterion  that  can  be 
applied  to  quantitatively  distinguish  between  isentropic 
CASC  results  is  required. 

In  the  absence  of  experimental  data,  it  is  common  analy¬ 
tical  design  practice  to  predict  aerodynamic  performance  by 
estimating  the  amount  of  diffusion  which  occurs  on  the  sur¬ 
face  of  a  given  airfoil  design.  Prince15,  McDonald17,  Szanca 
and  Schum3Q  address  the  topic  of  turbine  airfoil  diffusion 
and  concur  that  diffusion  control  is  a  significant  transonic 
blade  design  performance  criterion. 

Diffusion  or  flow  deceleration  causes  the  formation  of 
adverse  pressure  gradients  which  thickens  the  boundary  layer 
on  the  surface  of  the  blade.  Under  rapid  or  long  term  decel¬ 
erations,  the  flow  separates  completely.  The  possibility  of 
a  shock  wave  interaction  with  the  boundary  layer  of  a  dif¬ 
fusing  flow  will  almost  certainly  cause  a  flow  separation. 
Thus,  an  inviscid  prediction  of  the  amount  of  diffusion  which 
occurs  on  the  blade  suction  surface  is  an  indication  of  the 


tendency  for  the  flow  to  separate  and  thereby  significantly 
increase  pressure  losses. 


The  second  performance  criterion  that  will  be  applied 
to  the  various  blade  designs  is  summarized  in  the  following 
statement:  the  airfoil  shape  yielding  the  least  amount  of 
flow  overexpansion  and  subsequent  diffusion  is  the  better 
design.  This  flow  overexpansion  standard  is  applied  to  pro¬ 
posed  designs  through  use  of  the  following  equation; 

A  p  =  P  -  P  (13 

overexpansion  sexit  sminimum 


Evaluation  of  Existing  Design 

The  present  design  was  evaluated  in  terms  of  the  aero¬ 
dynamic  performance  criteria  developed  in  the  previous 
section  for  the  transonic  design  point.  Highlighting  the 
deficiencies  of  the  existing  airfoil  in  this  manner  demon¬ 
strates  the  need  for  design  of  a  transonically  well-behaved 
blade  shape . 

Since  there  is  no  transonic  experimental  data  for  the 
existing  turbine  airfoil,  the  CASC  flow  analysis  was  used 
to  obtain  predictions  of  the  flowfield  surrounding  the  exist¬ 
ing  airfoil  at  the  transonic  design  point.  In  this  case, 
the  CASC  input  $BB  namelist  is  exactly  the  same  as  those 
depicted  in  Appendix  A,  except  that  the  Table  V  thermodynamic 
data  are  substituted  in  lieu  of  the  verification  operating 
conditions.  For  choked,  transonic  flow  cases,  CASC  iterates 
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upon  the  upstream  inlet  Mach  number  (mass  flow  rate) ,  while 
holding  the  back  pressure  constant.  CASC  also  returns  exit 
static  pressure  for  either  blade  surface  that  satisfy  the 
isentropic  equations  and  boundary  conditions.  These  calcu¬ 
lated  exit  pressures  may  or  may  not  equal  each  other,  or  the 
back  pressure  depending  on  the  design. 

Figure  12  is  a  plot  of  the  pressure  distribution  re¬ 
sulting  from  the  evaluation  investigation.  Regions  of  rapid 
acceleration  can  be  observed  to  occur  on  both  the  pressure 
and  suction  surfaces.  Note  that  the  flow  over  the  suction 
surface  overexpands  to  a  very  high  Mach  number  (low  pressure) 
immediately  behind  the  sonic  throat  at  51  percent  of  axial 
chord.  After  this  initial  overexpansion,  the  velocity  (pres¬ 
sures)  on  the  suction  remain  relatively  constant.  Because 
the  flow  on  the  pressure  surface  does  not  exceed  a  Mach  num¬ 
ber  greater  than  one,  a  large  pressure  gradient  at  the 
trailing  edge  is  6.649  N/cmz  (9.63  psi) .  The  amount  of 
overexpansion  which  occurs  on  the  upper  surface  is  determined 
to  be  4.197  N/cmz  (6.08  psi). 

The  current  blade  behaves  poorly  because  too  much  turning 
occurs  on  the  suction  surface  after  the  throat  location. 

After  sonic  conditions  are  reached,  convex  surface  curvature 
causes  the  remaining  suction  surface  to  act  as  a  supersonic 
diverging  nozzle.  By  the  time  the  turning  becomes  small,  a 
high  velocity  has  been  attained  and  must  be  diffused,  sepa¬ 
rated  or  shocked-down  to  satisfy  the  exit  conditions  and 
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Figure  12  CASC  Prediction  of  the  Transonic  Design  Point 
CD  Distribution  for  the  Existing  Blade 
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trailing  edge  Kutta  boundary  requirements.  Finally,  lack  of 
a  convergent-divergent  flow  passage  prevents  the  flow  on  the 
pressure  surface  from  undergoing  supersonic  expansion  to  reach 
the  exit  static  pressure  level. 
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V.  Parametric  Design  Study 

The  geometry  of  the  existing  airfoil  typifies  most  sub¬ 
sonic  turbine  cascades.  A  turbine  airfoil  that  is  designed 
to  operate  at  subsonic  pressure  ratios  is  easily  identified 
by  the  following  geometric  characteristics;  a  convergent  flow 
passage,  significant  suction  surface  curvature  after  the  throat 
location,  a  maximum  blade  thickness  position  at  30-50  percent 
of  axial  chord,  and  a  moderate  stagger  angle.  The  net  effect 
of  these  blade  shape  trends  is  a  turbine  cascade  that  is  to¬ 
tally  unsuited  for  transonic  flow  environments. 

During  the  course  of  this  design  study,  the  observation 
was  made  that  a  well  designed  transonic  cascade  exhibits  two 
distinguishing  geometric  features.  First,  in  order  to  expand 
the  flow  on  the  pressure  surface  to  the  proper  exit  pressure 
level,  a  convergent-divergent  flow  passage  is  a  necessary  hall¬ 
mark  of  any  transonic  turbine  cascade.  Second,  practically 
all  of  the  suction  surface  flow  turning  must  occur  ahead  of 
the  throat  location  to  prevent  supersonic  overexpansion  and 
subsequent  diffusion  on  that  surface. 

Consequently,  the  twofold  objective  of  this  parametric 
design  study  is  to  identify  the  geometric  parameters  that 
have  a  major  influence  on  the  two  features  discussed  in  the 
previous  paragraph,  and  then,  to  determine  the  parameter 
combination  which  yields  the  maximum  amount  of  aerodynamic 
performance  at  the  design  point  operating  conditions. 
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Parametric  Design  Approacl 

Two  modes  of  data  input  to  the  BLDDES  module  were  dis¬ 
cussed  in  Chapter  II.  The  parametric  option  requires  utiliz¬ 
ation  of  a  BLDDES  subroutine  called  BSGEN,  which  analytically 
generates  blade  surface  coordinate  arrays  using  parametric 
input  in  lieu  of  digitized  data.  The  blade  section  generation 
(BSGEN)  calculation  procedure  involves  the  computation  of  a 
section  meanline  derived  from  a  fourth-order  polynomial  and 
a  thickness  distribution  derived  from  the  conformal  transform¬ 
ation  of  a  fifth-order  polynomial.  The  coefficients  of  these 
polynomials  are  governed  by  the  blade  shape  input  parameters. 

When  the  resulting  standardized  airfoil  definitions  are 
evaluated  by  the  aerodynamic  analysis  module,  CASC,  a  pressure 
distribution  prediction  for  the  original  set  of  blade  shape 
input  parameters  is  obtained.  Four  input  parameters  were 
found  to  have  a  substantial  effect  upon  the  airfoil  pressure 
distribution  predictions.  Thus,  the  design  study  was  limited 
to  variations  of  these  four  input  parameters.  The  influence 
of  each  input  parameter  upon  overall  aerodynamic  performance 
was  determined  by  generating  a  total  of  72  different  blade 
profiles,  and  then,  evaluating  the  associated  pressure  distri¬ 
bution  predictions. 

The  remaining  required  BLDDES  input  parameters  were  held 
constant  throughout  the  parametric  variation  study.  These 
fixed  parameters  generally  tended  to  relate  to  the  physical 
dimensions  of  a  given  airfoil  configuration  rather  than  the 
blade  profile  itself.  The  assigned  values  chosen  for  the 
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fixed  parameters  are  physically  compatible  with  the  static 
test  facility. 

Definition  of  Blade  Shape  Parameters 

The  primary  intent  of  this  section  is  to  identify  and 
adequately  describe  both  the  fixed  and  varied  blade  shape 
input  parameters.  A  complete  description  of  the  parametric 
input  variables  and  the  BSGEN  calculation  procedure  appears 
in  AFAPL-TR- 7  8 - 9  2  2  4  . 

Fixed  Parameters .  Table  VI  contains  a  listing  of  the 
parameters  that  were  assigned  fixed  values  during  the  para¬ 
metric  variation  study.  A  schematic  of  the  BSGEN  coordinate 
system  in  which  these  parameters  are  employed  is  presented 
in  Figure  13.  The  following  paragraphs  provide  a  brief 
description  of  the  parameters  tabulated  in  Table  VI. 

The  velocity  vector  diagram  information  that  is  consis¬ 
tent  with  the  test  facility  geometry,  is  represented  by  the 
first  two  Table  VI  entries.  BETA1  and  BETA2  are  the  specified 
gas  inlet  and  gas  exit  angles  respectively.  Together,  the 
two  quantities  fix  the  amount  of  flow  turning  to  be  accom¬ 
plished  in  the  cascade  passage. 

Throat  size  parameter,  CF,  influences  the  size  of  the 
flow  passage  geometric  throat  or  the  minimum  clearance  be¬ 
tween  blades.  The  throat  is  calculated  in  the  following 
manner; 

THROAT  =  CF  *  S  *  cos(BETA2)  (14) 
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Figure  13  Blade  Generator  Coordinate  System 


where  S  is  the  blade  spacing  (see  Figure  13) .  Initially, 
the  blade  exit  angle,  a2 ,  is  set  equal  to  the  gas  exit  angle, 
/?2*  If  the  throat  (minimum  blade  clearance)  calculated  by 
BSGEN  is  not  equal  to  the  desired  throat  represented  by 
equation  14,  a  new  blade  exit  angle  is  computed.  The  value 
chosen  for  CF  insures  that  the  mass  flow  rate  through  the 
test  facility  will  not  exceed  2.27  kgm/sec  (5  lbm/sec) . 

Although  the  parameter  C3  controls  the  intermediate 
curvature  of  the  section  meanline,  this  factor  has  a  minor 
effect  on  the  overall  blade  shape  when  compared  to  the  ef¬ 
fects  of  the  varied  parameters.  The  Table  VI  value  for  C3 
is  favorably  consistent  with  decreasing  the  suction  surface 
unguided  turning. 

The  angle  of  attack  to  the  blades  in  the  cascade  is 
specified  by  DELTAl.  Since  an  angle  of  attack  was  not  de¬ 
sired,  DELTA1  was  set  equal  to  zero.  Because  DELTAl  is  the 
angular  difference  between  the  tangent  line  to  the  section 
meanline  at  the  leading  edge  and  the  inlet  velocity  vector, 
the  following  relationship  exists  when  DELAT1  is  zero; 

dj  =  (15) 

where  c*i  and  @x  are  the  blade  inlet  and  gas  inlet  angles 
respectively. 

ISTK,  ITMODE,  and  IVANE  are  minor  blade  shape  indicators. 
These  quantities  were  defined  in  Chapter  III  and  the  fixed 
values  listed  in  Table  VI  are  consistent  with  the  previous 
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discussion.  NBLDS,  RLE,  RREF,  and  RTE  have  been  previously 
defined  in  Chapter  III.  The  number  of  blades  and  the  radii 
lengths  that  were  assigned,  define  a  cylindrical  blade  section 
with  a  cascade  blade  spacing,  S,  of  3.894  cm  (1.533  in).  The 
cylindrical  section  is  almost  rectalinear  because  of  the  large 
radial  values  chosen.  ZLE  and  ZTE  define  the  axial  width  of 
the  blades.  Note  that  the  blade  spacing  and  axial  width  for 
the  proposed  transonic  designs  is  the  same  as  that  for  the 
existing  cascade. 

Fifty  mils  is  the  minimum  tolerable  trailing  edge  thick¬ 
ness  required  for  advanced  turbine  blading  with  trailing  edge 
coolant  holes.  The  maximum  section  thickness  for  typical 
transonic  turbine  stators  is  usually  20  percent  of  the  axial 
chord  length.  This  general  rule-of-thumb  was  applied  to  fix 
TMAXX. 

Varied  Parameters .  The  four  blade  shape  input  parameters 
that  were  varied  are  listed  in  Table  VII.  The  input  parameter 
SF,  stagger  factor,  indirectly  controls  the  stagger  angle,  /?s, 
in  the  following  way; 

SF  =  tan  -  tan  0S  (16) 

tan  y3s  -  tan 

Hence,  the  relationship  between  stagger  factor  and  stagger 
angle  is  an  inverse  function,  i.e.,  increasing  SF  decreases 

The  stagger  angle  of  a  turbine  cascade  is  an  important 
transonic  design  parameter.  For  a  given  amount  of  flow 
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turning  (velocity  vector  diagram) ,  large  stagger  angles  tend 
to  reduce  surface  curvature  in  the  aft  portion  of  the  blades. 
Blades  with  high  stagger  angles  are  characterized  by  long, 
straight  tails  with  relatively  no  suction  surface  curvature 
after  the  throat  location.  The  primary  consequence  of  in¬ 
creasing  stagger  angle  is  to  shift  most  of  the  flow  turning 
to  the  foreward  portion  of  the  blades,  thereby,  enhancing 
transonic  performance.  Therefore,  the  stagger  factor  was 
treated  as  a  variable  design  parameter. 

The  input  parameter  Cl,  the  non-dimensional  axial  lo¬ 
cation  of  maximum  section  thickness,  was  also  treated  as  a 
design  parameter.  In  order  to  promote  higher  stagger  angles 
and  substantial  flow  turning  ahead  of  the  throat  location, 
most  of  the  blade  mass  (i.e.,  maximum  section  thickness), 
must  be  concentrated  towards  the  leading  edge  of  the  blade. 

The  maximum  section  thickness  location  along  with  stagger 
angle  primarily  determines  whether  or  not  the  flow  passage 
will  have  a  convergent-divergent  nature. 

Trailing  edge  taper  factor,  C2,  controls  the  trailing 
edge  included  angle,  6,  as  defined  by  the  following  equation; 

n  - 1  I  rrt/TMAXX-TE  \  q  I 

0  =  2  tan  VC2^ZTE---ZLEjCOS^  1  (17) 

The  trailing  edge  taper  factor  was  considered  to  be  a 
significant  design  parameter  because  trailing  edge  taper 
heavily  influences  the  throat  location  and  the  ratio  of  throat 
to  exit  area. 
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Finally,  the  leading  edge  bluntness  factor,  TI ,  is  re¬ 
lated  to  the  leading  edge  radius  of  curvature  as  shown  in 
Figure  13.  Bluntness  increases  with  increasing  TI  values. 

The  flowfield  behavior  over  the  midchord  portion  of  the  blades 
is  a  strong  function  of  the  leading  edge  bluntness.  Although 
the  parameter  TI  is  somewhat  dependent  on  the  maximum  section 
thickness  magnitude  and  location,  it  has  been  included  as  an 
independent  design  parameter. 

Transonic  Performance  Maximization 

Procedure .  The  four  blade  shape  input  parameters  (Cl, 

C2,  SF,  TI)  were  systematically  varied  to  generate  seventy- 
two  different  blade  design  versions  for  analysis  by  the  CAS C 
module.  Then,  the  resultant  pressure  distributions  were 
evaluated  in  terms  of  the  two  aerodynamic  performance  stan¬ 
dards  that  were  developed  in  Chapter  IV. 

The  specific  approach  that  was  used  to  maximize  the 
relationships  between  the  input  parameters  and  the  aerodyna¬ 
mic  performance  standards  was  a  two  phased  procedure.  Under 
phase  one,  the  input  parameter  combinations  which  yielded 
designs  with  trailing  edge  pressure  gradients,  off  design 
exit  static  pressures,  etc.,  were  identified  and  eliminated 
from  consideration.  This  objective  was  achieved  for  four 
separate  values  of  maximum  section  thickness  location  (Cl) , 
by  varying  the  other  three  input  parameters  until  a  unique 
design  was  found  that  satisfied  the  first  performance  cri¬ 
terion  at  the  design  point  exit  pressure.  The  second  phase 
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of  the  performance  maximization  process  was  concerned  with 
selecting  one  of  the  four  design  versions  from  phase  one. 

In  phase  two,  determination  of  the  superior  design  was  based 
upon  the  pressure  distribution  prediction  which  exhibited  the 
least  amount  of  suction  surface  diffusion. 

The  final  design  selection  is  presented  in  the  last  sec¬ 
tion  of  this  chapter.  Description  of  the  input  parameter/ 
aerodynamic  performance  relationships  that  were  discovered 
during  the  parametric  variation  study  is  the  subject  of  the 
next  two  subsections . 

Elimination  of  Trailing  Edge  Pressure  Gradients.  For  a 
given  value  of  maximum  section  thickness  location  (Cl) ,  the 
other  three  variable  blade  shape  input  parameters  (SF,  C2,  TI) 
were  individually  varied  to  determine  whether  or  not  any 
characteristic  effect  upon  pressure  distribution  shape  could 
be  attributed  to  a  specific  input  parameter.  Although  each 
of  the  three  input  parameters  did  not  exhibit  completely 
uncoupled  influences  upon  the  predicted  pressure  distributions, 
well ^defined  aerodynamic  trends  were  related  to  each  input 
parameter . 

The  stagger  factor,  SF,  primarily  determines  the  mag¬ 
nitude  of  the  pressure  gradient  across  the  blade's  trailing 
edge.  If  the  stagger  factor  is  too  large  (i.e.,  the  stagger 
angle  is  too  small)  for  a  given  set  of  operating  conditions, 
the  predicted  pressure  distribution  from  the  resulting  blade 
profile  will  feature  a  pressure  surface  that  is  underexpanded 
and  a  suction  surface  that  is  overexpanded.  Thus,  an  increase 
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in  the  stagger  factor  tends  to  decrease  the  convergent- 
divergent  nature  of  the  flow  passage.  A  limiting  case  is 
encountered  when  the  stagger  factor  becomes  great  enough  tc 
cause  the  flow  passage  to  become  entirely  convergent.  The 
opposite  extreme  condition  is  a  stagger  factor  that  is  too 
small,  which  yields  a  design  that  has  an  overexpanded  pressure 
surface  and  an  underexpanded  suction  surface.  Appendix  F 
contains  various  selected  design  versions  that  illustrate 
proper  and  improper  combinations  of  these  input  blade  shape 
parameters.  Design  version  two  is  an  example  of  a  blade 
contour  with  not  enough  stagger  angle,  while  design  version 
eighteen  typifies  the  case  with  too  much  stagger  angle. 

Once  the  appropriate  stagger  factor  is  found  such  that 
the  trailing  edge  pressure  gradient  is  eliminated,  the  input 
parameter  C2  is  varied  to  adjust  the  trailing  edge  taper 
ratio.  Although,  the  trailing  edge  pressure  gradient  may  be 
zero,  the  cascade  exit  static  pressure  to  which  both  surfaces 
are  expanded  may  not  necessarily  equal  the  design  back  pres¬ 
sure.  Design  version  twelve  (Appendix  F)  is  an  example  of  a 
design  case  where  the  trailing  edge  pressure  gradient  is  zero 
but  the  exit  static  pressure  is  considerably  lower  than  the 
design  back  pressure.  Decreasing  the  value  of  the  taper 
factor  has  the  effect  of  decreasing  the  trailing  edge  included 
angle  (see  Equation  17) .  The  ratio  of  the  flow  passage  throat 
to  exit  area  is  thereby  decreased  and  the  exit  static  pressure 
of  an  overexpanded  flow  is  raised.  Hence,  for  a  given  maxi¬ 
mum  thickness  location  and  stagger  factor,  proper  positioning 
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of  the  throat  location  with  reference  to  the  flow  passage 
exit  plane  is  accomplished  by  varying  the  taper  factor,  C2. 

The  last  variable  input  parameter  that  was  incorporated 
into  the  design  parameter  investigation,  is  the  leading  edge 
bluntness  factor,  TI.  The  combination  of  stagger  factor  and 
taper  factor  that  yields  a  design  with  zero  trailing  edge 
pressure  gradient  for  a  particular  exit  static  pressure  is 
still  not  a  unique  combination  until  the  parameter,  TI  is 
introduced.  The  characteristic  effect  on  predicted  pressure 
distribution  shape  that  is  associated  with  leading  edge 
bluntness  is  midchord  overexpansion  and  recompression  phe¬ 
nomena.  The  predicted  pressure  distribution  for  design 
version  fourteen  (Appendix  F)  shows  jagged  peaks  and  troughs 
that  are  commensurate  with  unstable  flow  accelerations  and 
decelerations  on  the  suction  surface.  A  larger  amount  of 
leading  edge  bluntness  alleviates  this  problem. 

Versions  fifty-two,  fifty-nine,  sixty-five,  and  seventy- 
two  are  the  design  versions  that  exhibited  zero  trailing  edge 
pressure  gradient  at  the  design  exit  static  pressure  with  no 
significant  midchord  overexpansions  and  recompressions.  These 
four  versions  respectively  correspond  to  decreasing  values  of 
maximum  section  thickness  location. 

Table  VIII  contains  a  complete  listing  of  the  variable 
input  parameter  combinations  used  to  generate  these  design 
versions.  Figures  14,  15  and  16  are  plots  of  SF,  C2  and  TI 
versus  the  maximum  section  thickness  location.  Cl,  for  the 
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Figure  14  Stagger  Factor  versus  Location 
of  Maximum  Section  Thickness 
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Figure  15  Taper  Factor  versus  Location 
of  Maximum  Section  Thickness 
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Figure  16  L.E.  Bluntness  Factor  versus  Location 
of  Maximum  Section  Thickness 


four  unique  design  versions  which  satisfied  the  first  aero¬ 
dynamic  performance  standard.  The  parametric  relationships 
expressed  in  these  figures  define  the  locus  of  designs  which 
yield  maximum  aerodynamic  performance  as  a  function  of  max¬ 
imum  section  thickness  location.  Version  seventy-two  was 
selected  as  the  best  design  and  therefore,  its  blade  contour, 
coordinates  and  predicted  pressure  distribution  appear  in 
the  last  section  of  this  chapter.  The  remaining  versions 
are  presented  in  Appendix  F. 

Minimizing  Suction  Surface  Diffusion .  The  design  ver¬ 
sions  listed  in  Table  VIII  are  the  superior  aerodynamic  blade 
contours  for  each  respective  maximum  section  thickness  lo¬ 
cation.  Selection  of  the  best  aerodynamic  profile  is  deter¬ 
mined  from  suction  surface  diffusion  considerations.  Dif¬ 
fusion  data  for  each  of  the  designs  is  tabulated  in  Table 
IX  and  plotted  in  Figure  17. 

Application  of  the  second  performance  crite"J.on  to  the 
diffusion  data  yields  version  seventy-two  as  the  best  overall 
aerodynamic  blade  design.  Although  the  Figure  17  diffusion 
trend  indicates  that  better  aerodynamic  designs  may  exist 
at  lower  maximum  section  thickness  locations,  design  versions 
with  maximum  section  thickness  location  ahead  of  the  20  per¬ 
cent  axial  chord  position  were  not  considered  due  to  practical 
mechanical  and  thermal  limits. 

Shifting  the  maximum  section  thickness  location  foreward 
while  maintaining  a  maximum  aerodynamic  performance  profile 
tends  to  increase  mechanical  and  thermal  stresses.  The 
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Figure  17  Normalized  Diffusion  versus  Location 
of  Maximum  Section  Thickness 
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effective  crossectional  area  for  interior  cooling  systems  and 
structural  integrity  continuously  reduces  with  decreasing  Cl 
values.  Not  only  does  the  overall  crossectional  area  become 
smaller,  but  the  area  distribution  within  the  profile  changes 
unfavorably,  i.e.,  the  area  in  the  front  portion  of  the  blade 
increases  while  the  area  in  the  already  thin  aft  portion 
decreases  even  further. 

Because  thermal  and  mechanical  analyses  were  not  in¬ 
cluded  within  the  scope  of  this  study,  an  arbritrary  thermal- 
mechanical  limit  was  imposed  upon  the  foreward  movement  of 
the  maximum  section  thickness  position.  The  Figure  14,  15, 
and  16  relationships  can  be  used  to  generate  new  designs 
should  a  thermal-mechanical  analysis  establish  that  this 
arbitrary  limit  is  not  sufficiently  restrictive. 

Final  Design 

On  the  basis  of  the  preceeding  analysis,  the  blade 
contour  which  yielded  the  maximum  amount  of  aerodynamic 
performance  is  design  version  seventy-two.  A  scaled  drawing 
of  the  blade  coordinates,  a  schematic  of  the  cascade  formed 
by  the  final  blade  profile,  and  the  predicted  pressure  dis¬ 
tribution  of  the  new  transonic  design  appears  in  Figures  18, 
19,  and  20.  The  blade  section  profile  coordinates  are  given 
in  Table  X.  These  coordinates  are  tabulated  with  respect  to 
an  origin  located  at  the  leading  edge  highlight  point,  and 
are  normalized  as  a  fraction  of  the  axial  width,  3.894  cm 
(1.533  in) . 
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Figure  19  Final  Transonic  Turbine  Cascade 
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Although  the  scaling  in  the  two  figures  is  different 
comparison  of  the  Figure  18  transonic  design  with  the  existing 
subsonic  profile  depicted  in  Figure  5  graphically  demonstrates 


the  dissimilarities  between  the  two  b]ade  contours.  The 
slightly  convergent-divergent  flow  passage  geometry  of  the  new 
transonic  turbine  cascade  is  evident  in  Figure  19.  Finally, 
the  inviscid  pressure  coefficient  distribution  prediction  of 
the  new  blade  profile  (Figure  20) ,  is  a  significant  improve¬ 
ment  over  the  current  design  (Figure  12) ,  in  terms  of  the 
Chapter  IV  performance  standards . 


VI .  Off-Design  Performance 

An  important  aspect  of  turbine  blade  design  is  the  pre¬ 
diction  of  cascade  performance  at  off-design  operating 
conditions.  Usually,  the  environment  under  which  the  turbine 
airfoil  actually  operates  is  significantly  different  from  the 
originally  intended  design  operating  conditions.  Therefore, 
the  new  transonic  turbine  stator  design  must  not  only  exhibit 
maximum  aerodynamic  performance  at  the  transonic  design  point 
specified  in  Chapter  IV,  but  it  should  exhibit  reasonable 
performance  over  a  wide  range  of  possible  test  environments. 

The  four  most  probable  off-design  scenarios  that  the 
new  design  is  anticipated  to  encounter  are  as  follows: 

1)  an  angle  of  attack  other  than  zero,  2)  a  different  gas 
total  temperature,  3)  cold  flow  tests,  and  4)  an  alternate 
pressure  ratio.  Hence,  the  predicted  performance  of  the 
proposed  airfoil  design  in  these  off-design  test  environments 
is  reported  in  the  following  sections  of  this  chapter.  In 
addition,  the  last  section  is  addressed  to  the  topic  of 
viscous  boundary  layer  influence  upon  aerodynamic  performance. 

Angle  of  Attack  Sensitivity 

The  final  turbine  stator  profile  presented  in  Chapter  V 
was  designed  for  an  angle  of  attack  equal  to  zero.  However, 
there  is  a  good  possibility  that  the  new  design  may  experience 
a  slight  angle  of  attack  if  the  gas  inlet  angle  in  the  test 
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facility  is  not  eight  degrees.  Although  the  Equation  8 
continuity  relationship  indicates  that  turbine  cascades  are 
insensitive  to  small  angles  of  attack,  the  effect  upon  per¬ 
formance  of  both  positive  and  negative  angles  of  attack  was 
analytically  investigated. 

Pressure  distribution  predictions  for  the  new  airfoil 
were  obtained  at  off-design  gas  inlet  angles,  ,  ranging 
from  zero  to  sixteen  degrees.  In  every  case,  the  results 
were  identical  to  the  pressure  coefficient  distribution 
plotted  in  Figure  20.  Thus,  the  recently  designed  transonic 
turbine  cascade,  sketched  in  Figure  19,  is  relatively  insen¬ 
sitive  to  small  angles  of  attack. 

Effect  of  Total  Temperature  Variations 

In  accordance  with  the  requirement  for  a  design  that  is 
based  on  an  extremely  high  gas  temperature,  a  design  point 
gas  total  temperature  of  3000°R  was  selected.  At  the  present 
time,  most  cooled  or  uncooled  turbine  stators  can  not  with¬ 
stand  these  temperatures  at  transonic  blade  loadings.  There¬ 
fore,  most  cascade  tests  in  the  immediate  future  will  probably 
occur  at  lower  gas  total  temperatures  until  superior  mater¬ 
ials  and/or  cooling  systems  are  developed. 

Thus,  the  effect  of  lowering  the  operating  gas  total 
temperature  to  values  of  2000°R  and  2500°R  v/as  investigated. 
Since  the  governing  continuity  and  momentum  equations  do  not 
have  explicit  temperature  terms,  the  effect  of  total  temper¬ 
ature  variation  was  expected  to  be  negligible.  This  outcome 
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was  realized  by  the  resulting  lower  temperature  pressure  dis¬ 
tributions  which  were  identical  to  the  Figure  20  pressure 
profile . 

Cold  Flow  Test  Conditions 

Cold  flow  tests  are  conducted  in  the  static  test  facility 
on  certain  occasions.  When  the  nature  of  the  experimental 
program  requires  use  of  low  temperature  instrumentation,  cold 
flow  testing  becomes  essential.  Therefore,  the  performance 
of  the  new  design  was  estimated  for  operation  in  this  type  of 
test  environment. 

Although  the  effect  of  total  temperature  variation  upon 
cascade  performance  was  reported  to  be  negligible  in  the 
previous  section,  cold  flow  test  conditions  marginally  effect 
aerodynamic  performance.  The  reason  for  this  apparent  dis¬ 
crepancy  is  due  to  the  fact  that  the  major  difference  between 
the  design  point  operating  condition  and  cold  flow  testing  is 
the  change  in  the  working  gas  medium,  not  the  total  tempera¬ 
ture  difference. 

During  the  design  point  operating  condition,  the  high 
gas  temperatures  are  generated  by  an  upstream  chemical  com¬ 
bustion  reaction  between  hydrogen  and  compressed  air.  When 
this  combustion  reaction  does  not  occur  (cold  flow) ,  the 
nature  of  the  working  gas  medium  considerably  changes.  Thus, 
it  is  the  specific  heat  ratio  and  molecular  weight  of  com¬ 
pressed  air  at  low  temperatures  which  causes  the  observable 
difference  in  aerodynamic  performance. 
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Figure  21  is  a  plot  of  the  pressure  distribution  predic¬ 
tion  for  a  cold  flow  test  condition  defined  by  the  following 
thermodynamic  quantities; 

TQ  =  300°K  ( 54  0°R) 

y  =  1.4 

R  =  286.9  - {1716.2  - ^o7t) 

sec2-°K  V  sec- -UR J 

The  isentropic  pressure  distribution  prediction  returned  by 
CASC  indicates  that  the  flow  on  the  suction  surface  will  tend 
to  undergo  an  increased  amount  of  diffusion  until  an  under¬ 
expanded  exit  pressure  of  13.563  N/cm2  (19.63  psia)  is 
reached.  Meanwhile,  the  flow  over  the  pressure  surface  is 
estimated  to  overexpand  to  an  exit  static  static  pressure 
of  12.405  N/cm2  (17.96  psia).  This  situation  does  not  occur 
in  an  actual  test  environment  because  loss  mechanisms, 
including  flow  separation  and  shock  wave  phenomenon,  equalize 
the  trailing  edge  pressure  gradient.  Hence,  overall  aero- 
cynamic  performance  degrades  during  cold  flow  operation.  On 
the  basis  of  the  design  parameter  relationships  discussed  in 
Chapter  V,  this  performance  problem  can  be  partially  alleviated 
by  slightly  decreasing  the  cascade  stagger  angle,  i.e., 
increase  the  stagger  factor,  SF. 

Alternate  Pressure  Ratio 

According  to  the  qualitative  functional  relationship 
expressed  by  Equation  9,  the  pressure  coefficient  distri¬ 
bution  is  heavily  dependent  upon  the  cascade  pressure  ratio. 
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Figure  21  Off-Design  CD  Distribution  Prediction 
(Cold  Flow  Test  Condition) 
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Therefore,  both  subsonic  and  transonic  off-design  pressure 
ratios  were  analytically  applied  to  the  new  design.  In  terms 
of  significant  loss  in  cascade  performance,  the  results  of  this 
off-design  pressure  ratio  investigation  are  the  most  dramatic. 

Subsonic .  An  evaluation  of  the  existing  design  deter¬ 
mined  that  subsonic  blade  geometries  exhibit  poor  aerodynamic 
performance  at  transonic  flow  conditions.  In  Chapter  V,  the 
geometric  qualities  which  distinguish  a  subsonic  turbine  shape 
from  a  transonic  contour  were  discussed.  This  section  is 
concerned  with  presenting  the  performance  characteristics  of 
the  new  transonic  design  with  respect  to  a  subsonic  pressure 
ratio. 

When  the  recently  developed  transonic  design  is  analy¬ 
tically  subjected  to  a  subsonic  pressure  ratio,  the  isentropic 
pressure  coefficient  distribution  displayed  in  Figure  22  is 
generated.  Large  amounts  of  flow  overexpansion  are  predicted 
on  both  the  pressure  and  suction  surfaces.  In  fact,  the  ex¬ 
istence  of  a  supersonic  bubble  on  the  suction  surface  is 
estimated  to  occur  at  30-50  percent  of  axial  chord.  Along 
the  divergent  portion  of  the  flow  passage,  flow  on  both  sur¬ 
faces  diffuses  drastically.  Again,  the  situation  depicted 
in  Figure  22  is  not  realistic  at  the  trailing  edge  because 
flow  separation  and  other  loss  mechanisms  will  cause  the 
flows  over  the  aft  portion  of  both  surfaces  to  separate, 
rather  than  isentropically  diffuse  to  meet  exit  static 
pressure  boundary  conditions. 
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Figure  22  Off-Design  Cp  Distribution  Prediction 
(Subsonic  Pressure  Ratio) 
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Transonic.  Figures  23  and  24  show  the  effect  of  re¬ 
spectively  increasing  and  decreasing  the  pressure  ratio  while 
maintaining  a  transonic  operating  condition.  The  cascade 
pressure  ratio  can  be  increased  by  either  increasing  the 
total  pressure,  PQ/  or  decreasing  the  back  pressure,  Pb*  The 
opposite  relationships  apply  for  decreasing  the  pressure 
ratio. 

Increasing  the  cascade  pressure  ratio  above  the  design 
point  enhances  the  tendency  of  a  negative  trailing  edge 
pressure  gradient  to  be  predicted,  i.e.  ,  the  suction  surface 
overexpanded  and  the  pressure  surface  underexpanded.  A 
decrease  in  the  cascade  pressure  ratio  below  the  design 
point  has  the  opposite  effect;  positive  trailing  edge  pres¬ 
sure  gradients  are  predicted.  These  tendencies  are  demon¬ 
strated  in  Figures  23  and  24  respectively.  Stagger  angle 
can  be  adjusted  to  correct  these  adverse  aerodynamic  effects 
in  the  following  manner;  reduce  stagger  angle  to  compensate 
for  negative  pressure  gradients,  increase  stagger  angle  to 
compensate  for  positive  pressure  gradients. 

Effect  of  Viscous  Boundary  Layer 

The  final  transonic  design  is  based  upon  maximizar.ion 
of  aerodynamic  performance.  However,  the  flowfield  calcu¬ 
lation  procedure  and  thermodynamic  standards  utilized  to 
achieve  this  result  are  predicted  upon  isentropic,  inviscid 
theory.  In  an  actual  flow  situation,  viscous  boundary  layers 
are  present.  The  overall  effect  of  boundary  layer  development 
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Figure  23  Off-Design  Cp  Distribution  Prediction 

(Transonic  Pressure  Ratio  above  Design  Point) 
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Figure  24  Off-Design  Cp  Distribution  Prediction 

(Transonic  Pressure  Ratio  below  Design  Point) 
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is  the  existence  of  an  additional  displacement  thickness. 
Since  the  flow  area  is  reduced,  higher  velocities  result  for 
a  given  mass  flow  rate. 

This  situation  was  simulated  by  considering  a  design 
version  with  a  slightly  larger  thickness  distribution  along 
the  latter  portion  of  the  blade.  Although  the  exit  static 
pressures  at  the  trailing  edge  remain  equal,  the  overall  exit 
static  pressure  decreases  below  the  design  back  pressure  by 
2  percent.  Thus,  boundary  layer  induced  effects  upon  perfor¬ 
mance  are  relatively  insignificant. 
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VII.  Conclusions  and  Recommendations 

A  turbine  stator  profile  that  is  analytically  predicted 
to  exhibit  good  aerodynamic  performance  over  a  wide  range 
of  transonic  operating  conditions  has  been  designed  fcr  the 
AFAPL  Heat  Transfer  Test  Facility.  The  design  of  the  new 
turbine  airfoil  section  was  based  entirely  upon  analytical 
methods.  THis  is  a  revolutionary  approach  to  the  turbine 
blade  design  problem.  During  the  course  of  this  design 
study,  a  number  of  significant  conclusions  were  formulated. 

First,  the  analytical  TDS  computer  codes  that  were  used 
to  generate  blade  profiles  and  perform  cascade  flow  calcu¬ 
lations  were  verified  for  accuracy  and  consistency.  Analytical- 
experimental  correlation  of  TDS  pressure  distribution  predic¬ 
tions  with  existing  test  facility  cascade  data  demonstrated 
the  ability  of  the  computer  programs  to  compute  heat  transfer 
facility  cascade  flowfields  successfully. 

A  set  of  aerodynamic  performance  standards  to  be  applied 
to  isentropic  cascade  flowfield  solutions  has  been  developed. 
Although  experimental  performance  data  is  conclusive,  appli¬ 
cation  of  isentropic  standards  to  state-of-the-art  analytical 
pressure  distribution  predictions  is  less  time  consuming, 
more  cost-effective,  and  reasonably  valid. 

A  well  designed  transonic  turbine  stator  must  have  a 
convergent-divergent  flow  passage  with  most  of  the  flow 


turning  occuring  ahead  of  the  throat  location.  Important 
geometric  parameters  for  the  design  of  transonic  turbine 
airfoils  have  been  identified.  Transonic  turbine  stators 
require  large  stagger  angles  compared  with  typical  subsonic 
designs.  The  stagger  angle  is  the  critical  geometric  para¬ 
meter  that  determines  whether  or  not  the  flow  is  equally 
expanded  on  both  the  pressure  and  suction  surfaces.  In  ad¬ 
dition,  the  location  of  maximum  section  thickness  should  be 
positioned  ahead  of  the  25  percent  axial  chord  location  in 
order  to  minimize  diffusion  tendencies. 

The  results  of  the  off-design  investigation  indicate 
that  the  performance  of  the  final  transonic  design  is  not 
seriously  compromised  at  ultimate  angles  of  attack,  and 
total  temperature  values.  However,  small  changes  in  the 
operating  pressure  ratio  will  significantly  change  the 
pressure  distribution  of  the  new  design.  From  an  aerodyna¬ 
mic  performance  or  diffusion  control  viewpoint,  operation 
of  the  final  transonic  design  at  subsonic  pressure  ratios 
would  be  intolerable. 

The  scope  of  this  transonic  turbine  stator  design  study 
was  limited  to  maximization  of  aerodynamic  performance. 
Further  study  is  needed  to  develop  a  thermal  and  mechanical 
analysis  suitable  for  incorporation  into  the  performance 
maximization  process. 

A  second  area  requiring  further  study  is  development 
of  an  improved  aerodynamic  analysis  which  incorporates  non- 
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isentropic  flow  phenomena.  For  instance,  formulation  of 
viscous  trailing  edge  wake  model  might  significantly  improve 
the  predicted  flowfield  calculations. 

The  final  recommendation  is  that  the  turbine  stator 
design  presented  in  this  report  be  built  and  tested  in  the 
Air  Force  Aero  Propulsion  Laboratory  Heat  Transfer  Test 
Facility.  The  results  of  these  experimental  tests  could 
provide  valuable  data  regarding  the  analytical  concepts 
used  to  generate  this  design. 
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APPENDIX  A 

CASC  Input  Files  for 
Test  Cases  One  and  Two 
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APPENDIX  B 

Derivation  of  the  Non-Dimensional 
Turbine  Cascade  Continuity  Equation  (Eq.  8) 


1 
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Equation  B-l  is  the  simple  continuity  balance  which 
states  that  the  mass  flow  through  the  inlet  plane  is  equal 
to  the  mass  flow  out  the  exit  plane.  This  is  equation  seven 
in  the  text  (Chapter  III) . 


Aj  =  m2  (B-l) 

The  subscripts  1  and  2  refer  to  the  turbine  cascade 
inlet  and  exit  stations  respectively  (see  Figure  6) .  Equa¬ 
tion  B-l  can  be  expressed  in  the  Figure  6  terminology; 

PlAjVzj  =  P2A2VZ2  (B-2 ) 


where 

p  =  local  density 
A  =  flow  area 
Vz  =  axial  velocity 


It  should  be  readily  apparent  that  the  following  def¬ 
initions  are  valid  for  a  turbine  cascade  geometry  (see  Figure 
6  control  volume) : 


A  2  =  S*L 

(B-3) 

(cos  /?!)  V 1 

(B-4 ) 

(cos  p2)  v2 

(B-5) 

where 

S  =  cascade  pitch  (blade  spacing) 
L  =  characteristic  unit  length 
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f$l  =  gas  incidence  angle 
f$2  =  gas  exit  angle 
Vi  =  inlet  velocity 
V2  exit  velocity 

Substituting  Equations  B-3-  B-4,  and  B-5  into  Equation 
B-2  and  dividing  through  by  the  equal  area  term  yields  the 
following  mass  flux  expression; 

Pi  (cos  fii)Vx  =  p2(cos  y92)v2  (B-6) 


From  gas  dynamic  theory,  velocity,  density  and  tempera¬ 
ture  can  be  represented  in  the  following  manner: 


V  =  MyJyRT 

(B-7) 

1 

p  v-i  i ~y 

P  ’  R?0(1  +  2  M2» 

(B-8) 

T  =  T0/(l  +  —■  M2) 

(B-9) 

Substitution  of  Equation  B-9  into  B-7  yields: 


(B-10) 


After  the  appropriate  subscripts  are  introduced,  Equations 
B-10  and  B-8  can  be  inserted  into  Equation  B-6  to  produce  the 
following : 
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APPENDIX  C 

CASC  Output  Files  for 
Test  Cases  One  and  Two 
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APPENDIX  D 

Pressure  Coefficient  Distributions 
for  Various  Combinations  of  /?  j  and  Mx 
(First  Test  Case) 
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FRACTN  AXIAL  CHORD 


Figure  25  Case  One  Cp  Distribution  Correlation 
(Mi  =  0.2775) 
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FRflCTN  AXIAL  CHORD 


Figure  26  Case  One  Cp  Distribution  Correlation 
(Mx  =  0.2813) 
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FRACTN  AXIAL  CHORD 


Figure  27  Case  One  Cp  Distribution  Correlation 
(Mj  =  0.2850) 
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FRflCTN  RXIRL  CHORD 


Figure  28  Case  One  Cp  Distribution  Correlation 
(M!  =  0.2900) 
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Figure  29  Case  One  C_  Distribution  Correlation 
(Mi  =  0.2950) 


120 


APPENDIX  E 

Pressure  Coefficient  Distributions 
for  Various  Combinations  of  /S  x  and  Mj 
(Second  Test  Case) 
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Figure  30  Case  Two  Cn  Distribution  Correlation 
(Mi  =  0.2950) 
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PS/PT0TRL 


Figure  31  Case  Two  Cp  Distribution  Correlation 
(Mi  =  0.2970) 
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FRflCTN  flXIRL  CHORD 


Figure  32  Case  Two  Cp  Distribution  Correlation 
(Mx  =  0.2980) 
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Figure  33  Case  Two  Cp  Distribution  Correlation 
(Mi  =  0.3000) 
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Figure  34  Case  Two  Cp  Distribution  Correlation 
=  o.30in> 
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APPENDIX  F 

Selected  Design  Versions 
and 

Associated  Cp  Distribution  Predictions 
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Figure  35  Design  Version  Two 
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Figure  37  Design  Version  Twelve 
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Figure  38  Design  Version  Twelve  Predicted  Cp  Distribution 
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Figure  39  Design  Version  Fourteen 
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Figure  41  Design  Version  Eighteen 
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Figure  43  Design  Version  Fifty-Two 
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Figure  44  Design  Version  Fifty-Two  Predicted  Cp  Distribution 
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Figure  45  Design  Version  Fifty-Nine 
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Figure  46  Design  Version  Fifty-Nine  Predicted  Cp  Distribution 
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Figure  47  Design  Version  Sixty-Five 
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